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Abstract 

We prove a conjecture made in |C1| : given an n-body central configu- 
ration Xo in the euclidean space E of dimension 2p, let ImJ' be the set of 
ordered real p-tuples {vi, vi^ - ■ ■ , i^p} such that {iiui, ±ii/2, • • • , ±ii/p} is 
the spectrum of the angular momentum of some (periodic) relative equi- 
librium motion of Xq in E. Then ImT is a convex polytope. The proof 
consists in showing that there exist two, generically {p — l)-dimensionaI, 
convex polytopes Vi and in W such that Vi C ImJ^ C 1^2 and that 
these two polytopes coincide. 

Vi, introduced in [CI] , is the set of spectra corresponding to the her- 
mitian structures J on E which are "adapted" to the symmetries of the 
inertia matrix So; it is associated with Horn's problem for the sum of 
px p real symmetric matrices with spectra ct_ and a+ whose union is the 
spectrum of So. 

7^2 is the orthogonal projection onto the set of "hermitian spectra" of 
the polytope V associated with Horn's problem for the sum of 2p x 2p real 
symmetric matrices having each the same spectrum as So- 

The equality Vi = V2 follows directly from a deep combinatorial 
lemma, proved in [FFLPj . which characterizes those of the sums C = A+B 
of two 2p X 2p real symmetric matrices A and B with the same spectrum, 
which are hermitian for some hermitian structure. 

1 Origin of the problem: A^^-body relative 
equilibria and their angular momenta 

We recall here the results of [AC I ICll IC2] which are needed in order to 
understand the mechanical origin of the purely algebraic conjecture solved 
in the present paper: given a configuration xq = (ri,--- ,rjv) € E'^ of 
A'^ punctual positive masses in the euclidean space E, a rigid motion of 
the configuration under Newton's attraction is a motion in which the mu- 
tual distances \\fi — rj\\ between the bodies stay constant. It is proved 



in [AC] (see also [C2| ) that such a motion is necessarily a relative equi- 
librium. This implies that the motion takes place in a space of even 
dimension 2p, which can be supposed to coincide with E, and that, in 
a galilean frame fixing the center of mass at the origin, it is of the form 
x{t) = (e^'^fi, e^*f2, ■ ■ ■ , e^*riv), where n is a 2p x 2p-antisymmetric en- 
domorphism of the euclidean space E which is non degenerate. Choosing 
an orthonormal basis where il is normalized, this amounts to saying that 
there exists a hermitian structure on the space E of motion and an or- 
thogonal decomposition E = C' = C''^ x ■ • • x C*'' such that 

x(t) = (xi(t), - ■ ■ , Xr{t)) = (e'"^*!!, ■ • ■ , e'"''*Xr), 

where Xm is the orthogonal projection on C*™ of the A^-body configuration 
X and the action of e'"""' on Xm is the diagonal action on each body of 
the projected configuration. Such quasi-periodic motions exist only for 
very special configurations, called balanced configurations (see [ACl IC2| 
for their characterization). The most degenerate balanced configurations 
are the central configurations for which all the frequencies LUi are the same; 
this means that 57 = ljJ, with J a hermitian structure on E, and the 
motion is 

x{t) = (flit), - ■ ■ ,rjv(t)) = e"^XQ = (e'"Vi, ■ • • , e"^V]v) 

in the hermitian space E = C; in particular, it is periodic. In a space of 
dimension at most 3, E is necessarily of dimension 2 and the configuration 
of any relative equilibrium is central. 

Given a configuration x = (ri,--- ,fjv) and a configuration of ve- 
locities y = X — {vi,--- ,vn), both with center of mass at the origin: 
XlfcLi "T-fc'^fe — X^fcLi ^fc'^'fc — 0, the angular momentum of {x,y) is the 
bivector C = X^^^ rntrk /\ Vk- If we represent x and y by the 2p x N 
matrices X and Y whose ith column are respectively made of the com- 
ponents (rii, • • • , r2pi) and {vu, ■ • ■ , V2pi) of ri and Vi in an orthonormal 
basis of E and if A'l — diag(mi,--- ,mN), this bivector is represented 
by the antisymmetric matrix (we use the french convention *X for the 
transposed of X) 

N 

C — —XM^Y -\- YM^X with coefficients Cij — "^^mki—rikVjk + rjkVik)- 

k = l 

The dynamics of a solid body is determined by its inertia tensor (with 
respect to its center of mass), represented in the case of a point masses 
configuration X by the symmetric matrix 

N 

S — XM^X with coefficients Sij = mkrjkrjk, 

k = l 

whose trace is the moment of inertta of the configuration x with respect 
to its center of mass. In particular, the angular momentum of a relative 
equilibrium solution X{t) = e*^A"o is represented by the antisymmetric 
matrix C = Sofi -I- fiSo, where 5*0 = XqM^Xo- Restricting to the case of 
central configurations, that is f7 = ojj, and making a; = 1, we consider 
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in what follows the spectrum of J-skew-hermitian matrices of the form 
SoJ + J So or, what amounts to the same, the spectrum of J-hermitian 
matrice^of the form J^^SoJ + So- 

In the following, we identify E with R^*' by the choice of some orthonormal 
basis. R^*" is endowed with its canonical basis a = (0, • • ■ , 1, ■ • • ,0) and 
its canonical euclidean scalar product x ■ y = X^Si ^JiJ/i,' this allows iden- 
tifying linear endomorphisms of E — R^'' and 2p x 2p matrices with real 
coefficients. When we say that J is a hermitian structure, we mean that 
the standard euclidean structure is given and that J is a complex structure 
which is orthogonal. 



2 The frequency map 



We recall the definition, given in [CTj . of the frequency map T from the set 
of hermitian structures on R'^'' to the positive Weyl chamber C R*": 
given some 2p x 2p real symmetric matrix So, we consider the mapping 
J I— > J~^SoJ + So from the space of hermitian structures on E to the 
set of 2p X 2p real symmetric matrices. We are only interested in the 
spectra of these matrices, hence choosing an orientation for J is harmless 
and we shall consider only those of the form J — JqR, where Jo is 

the "standard" structure defined by Jq = ( ^ ^\ and R G SO{2p). 



Id 

Two elements R' and R" of SO{2p) defining the same J if and only if 
there exists an element U £ U{p) such that R!' = UR! , it follows that 
the space of (oriented) hermitian structures is identified to the homo- 
geneous space U{p)\SO{2p). The symmetric matrix J~^SoJ + So is J- 
hermitian, that is, it commutes with J. This implies that its spectrum 
is real, of the form {i/i, !^2, • ■ • ,i^p} if considered as a p x p complex ma- 
trix (for the identification of R^'' to C defined by J) and of the form 
{{t^i, i'2, - ■ ■ , fp), {i^i, V2, - ■ ■ , i^p)} if considered as a 2p x 2p real matrix 
(see the next section for the trivial proof). 

Definition 1 The frequency mapping 

: U{j>)\SO{2p) -> = {(!/i, ■ ■ ■ i^p) G R^ 1^1 > ■ • ■ > fp} 

associates to each hermitian structure J the ordered spectrum {v\, • • • , Vp) 
of the J-hermitian matrix J~^SoJ + So. 



3 Hermitian spectra 

Lemma 1 Let C : M.'^^ — >■ R^'' be a symmetric endomorphism. The fol- 
lowing assertions are equivalent: 

1) There exists a hermitian structure J = R~^ JqR such that C be 
,J -hermitian (i.e. JC = CJ); 

2) The spectrum cy(C) of C is of the form 

o"(C) = {{'^i,'^2,- ■ ■ ,fp),(t'i,i^2,--- ,fp)}- 

^Notice that this is the same as the spectrum of the Jo-hermitian matrix S = J^^ SJo -\-S, 
where S = RSoR~^, which was considered in [CI) . 
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Proof. Let J = RT^JqR; the mapping C is J-hermitian if and only if 
RCR~^ is Jo-hermitian. This is equivalent to the existence of an isomor- 
phism U € U(p) C SO{2p) such that 

U~^RCR~^U = diag ((;/!, •• • ,i>p),{i/i,- ■ ■ ,i^p)). 

Conversely, the identity 

R~^CR = diag((i^i, • • • ,i'p), {vi,-- • jt'p)) 

implies the commutation of R~^CR with Jo and hence the commutation 

of C with J = 7?- Voi?. 

Notations. We shall call hermitian the spectra of this form and the 
diagonal the linear subspace A of VFj^ defined by 

A = {(/Xl > M2 > • • • > M2p), Ml = A*2, At3 = M4, ■ ■ ■ , /"2p-l = M2p}- 

Hence the ordered hermitian spectra are the ones belonging to A. 

4 Two convex polytopes 

Let So : R^*" — ^ R^*" be a symmetric endomorphism with spectrum 

a{So) = {fTl > 0-2 > • • • > CT2p}. 

To So we associate the subsets Vi and V2 of R^ (in fact of the positive 
Weyl chamber Wj+ of ordered real p-tuples), defined as follows: 

1) Vi is the set of ordered spectra 

a-(c) = {i/i > f2 > • • • > Vp} 

of symmetric endomorphisms c of of the form c = a + b, where a : 
M'' — ^ M.'' and 6 : M'' — ^ M'' are arbitrary symmetric endomorphisms with 
respective spectra 

cr(o) = CT- := {(Ti, era, • • • , cr2p-i}, cr(6) = cr+ := {(72, C4, • • • , cr2p}; 

2) ■P2 is the set of p-tuples {vi > 1/2 > ■ ■ ■ > Vp} such that 

{(z/i, 1/2 • • • ,Up),{ui,i'2--- , i-p)} 

is the spectrum of some symmetric endomorphism C of R^** of the form 
C = A + B, where A : R'^p ^ R^p and B : R^" R^** are arbitrary 
symmetric endomorphisms with the same spectrum 

= a(S„). 

In other words, identifying canonically the diagonal A with R^, one can 
write 

■P2 = P n A, 

where V is the (2p — l)-dimensional Horn polytope which describes the 
ordered spectra of sums C = A + B of 2p x 2p real symmetric matrices 
A, B with the same spectrum as So- 
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Lemma 2 Vi and are both contained in the hyperplane of W with 
equation 

i=l j = l 

Moreover, Vi andV2 are both {jp — 1)- dimensional convex polytopes and 

Vi C ImT C ■p2. 

Proof. The first identity comes from the additivity of the trace function. 
The fact that both Vi, V and hence 1^2 = V C] A, are convex polytopes 
is a general fact coming from the interpretation of the Horn problem as 
a moment map problem. Finally, the second inclusion comes from the 
very definition of and the first comes from Lemma 1 and the following 
identity, where (t_ and a+ axe considered as p x p diagonal matrices and 
P G SO(p): 

a+) \p J \0 a+) \p J 

CT- + p^^CT+p \ 

pa^p^^ + o-+y 

Remark. The choice of the partition cr = cr_ U cr+ of cr is dictated by 
the following theorem, which proves that any other partition of cr into two 
subsets with p elements will lead to a smaller polytope Vi : 

Theorem 3 C |FFLP| Proposition 2.2) Let A and B be p x p Hermi- 
tian matrices. Let ai > a"2 > • • ■ > (J2p be the eigenvalues of A and B 
arranged in descending order. Then there exist Hermitian matrices A and 
B with eigenvalues ai > erg > ■ • ■ > a2p-i and o"2 > (T4 > • • • > (J2p 
respectively, such that A + B = A + B . 

This was used in [CI] to prove that Vi — ImA is the image under the 
frequency map J- of the adapted hermitian structures, i.e. those J which 
send some p-dimensional subspace of M.^'' generated by eigenvectors of So 
onto its orthogonal. 

5 The projection property 

In this section, we prove the 

Theorem 4 The two polytopes Vi and 1^2 coincide. 

Corollary 5 ImT — Vi — ImA. In other words, the image by the 
frequency map J- of the adapted structures is already the full image IraJ- . 

We need recall the inductive definition of the Horn inequalities which 
define the Horn polytopes (see [F]). For a sum a + 6 = c of symmetric 
p X p matrices with respective (ordered in decreasing order) spectra 

a = (qi, • • • ,Qp), /3 = (^1, • • • ,^p), 7 = (71, • • • ,7p). 
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they read 



CIJK) Vr < p, V(/, J, K) e T^, J2^^^J2'^^+Y1 

where (notation of |F], noted LR^ by reference to Littlewood-Richardson 
coefficients in |FFLP| ') is defined as foUows: let Ur be the set of triples 
{I, J, K) of subsets of cardinal r of {1, 2, • ■ • ,p} such that 



ie/ j6J fcGA' 



r{r- + 1) 



Then set Tf = Uf and define recursively by 

rpp 



(/, J, K) G U?,-is < r, V(F, G, i?) G TJ, 

An immediate computation gives the following 
Lemma 6 Let 



h 


= (2ii 


- I,2i2 - 1,- 




l,2ji,2j2, 


•• ,2jV), 




= (2ii 


- 1,2*2 - 1,- 


■ ■ , 2ir - 


l,2ji,2j2, 


■■ ,2jr), 


K2 


= (2fci 


- I,2fci,2fc2 


- I,2fc2, 







Then {I2,J2,K2) G U^^ . 
Proof. It suffices to check that 



^(2i-l) + ^(2j) 



= ^(2fc-l)+^2fc + 



2r(2r + 1) 



Now, comes the key fact: 

Theorem 7 ([FFLP], lemma 1.18) For any triple (/, J, K) in , the 

triple {I2, J2, -ft'2) is in T^^ 

The proof of this theorem, which concerns the so-called "domino-decompo- 
sable Young diagrams" , is based on a version of the Littlewood-Richardson 
rule due to Carre and Leclerc ICLI. 



It implies that, for any a sum A + B — C of real symmetric 2p x 2p 
matrices with respective (ordered in decreasing order) spectra 

a = {&!,■ ■ ■ ,&2p), /3 = (/3i, ■ • ■ ,/32p), 7 = (71, • ■ • ,72p), 

(* I2,J2, K2) holds, that is 

E (72fe-l + 72fc) < E('^2»-l + P2i-l) + E('^2j + /32j). 

fegx 16/ jeJ 

In particular, if 

we get that 



E 



a = P = O = ((Jl, (72 
72fc-l + 72fc 



, cr2pj 



< E'^2i-i + E°"2i- 
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Note that the mapping 



(71,72,- •■ ,72p-l,72p) ( 



71 + 72 71 + 72 
2 ' 2 



72p-l + 72p 72p-l + 72p 

2 ' 2 



) 



is the orthogonal projection of the ordered set (71,72, • • • ,72p-i,72p) on 
the diagonal A of K^'' defined by the equations 71 = 72, • • • , 72p-i = 72p, 
that is on the subset of "hermitian" spectra. Hence a paraphrase of the 
above theorem is 

Theorem 8 Let C — A + B be the sum of two 2p x 2p real symmetric 
matrices with the same spectrum {ai > > ■ ■ ■ > o"2p}- 
If {1^1 > • • • > I'p} *s the orthogonal projection on the diagonal A = of 
the spectrum {71 > 72 > ■ • ■ > 72p} of C , that ts if Uk — llh^z^tJllL^ ffig 
triple of ordered spectra 

a = ((Jl, as, • ■ • , (T2p-l), P = (0-2,0-4, • • • ,02p), 7 = (i^i, i/2, • • • , fp) 
satisfies the Horn inequality (*I,J,K). 

This implies the follo-wing extremal property of the subset of "hermitian" 
spectra: 

Corollary 9 The orthogonal projection on the diagonal A of the (2p— 1)- 
dimensional Horn polytope V C M^'' associated with the spectra 



coincides with the {p — 1)- dimensional Horn polytope Vi G W associated 
with the spectra 



In particular, the intersection P2 = "P n A corresponding to the hermitian 
spectra, that is those such that 71 = 72, - - - ,72p-i = 72p, coincides -with 
Vi- Indeed, 1^2 coincides with the projection of V, which itself coincides 
with Vi- 

Remark. The equality ImT = 1^2 implies the following 

Corollary 10 Let C : R^^ -> R^*" be the sum C ^ A+B of two symmetric 
endomorphisms A and B with the same spectrum o-{A) — (t{B) = a{So). 
Then C is J-hermitian for some hermitian structure J on R^'' if and 
only if it is conjuguate by an element of SO{2p) to a matrix of the form 
So + J~^SoJ, where J is a hermitian structure on H?^ . 

Note. The proof of the above results has been written by the first author 
after he was convinced by the numerical experiments made by the second 
author that the equality Vi = 7^2 was plausible when p = 3 and more 
precisely that not only the intersection 1^2 = Pn A but also the orthogonal 
projection of the Horn polytope P on A was contained in Vi after the 
canonical identification of A with R*". This led first to a proof when p = 2 
or 3, obtained by coping directly with Horn's inequalities and then to the 
discovery that the general case followed from a lemma which turned out 
to be exactly the lemma 1.18 of [FFLP] . The numerical experiments are 
described in the next section. 



a{A) = a{B) = {01 > cr2 > - ■ - > o-2p} 



a{a) = (cti, 0-3, • • ■ , o-2p_i), o(fe) = (02, 04, • • • , 02p). 
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6 Numerical experiments 



The numerical checking of the conjecture that JmJ-" = Vi, was made 
on the matrix So = ^diag {13, 8, 5, 3, 2, 1}, whose spectrum satisfies the 
inequalities in [CI) (section 8). We wrote a program in TRIP [GLllj 
producing different rotation matrices R £ SO{2p) in a random way. 
Starting from the canonical basis ^ = {^i,...,^m}, rn = p{2p — 1), 
of so(2p), we created a list containing the m one-parameter subgroups 
Gi{t) = e*^' C SO{2p). We created a second list of m random values 
and a random permutation [1, 2, . . . , m] — >• [ii, 12, ... , im]- The ran- 
dom rotation matrix was defined as 

m- 

R=Y[G,^{tj), 0<U<2n. (1) 
3=1 

The program which plots Vi is very simple (the fact that we replaced the 
conjugation of Jq by the conjugation of is immaterial and comes from 
the formulation of the conjecture in [Cl| ): 



create So and Jo 
for z = 1 to Nmax do 

create R and R~^; 

S = RSqR ^ ] 

C = S — JqSJq', 

1st = eigenvalues (C) ; 

plot ( 1st [5], 1st [3], lst[l] ); 
end for. 

We have assigned the value Nmax = 25000 obtaining the results shown 
in Figure [T] The figure shows also the simplex 71 -|- 72 + 73 = 1 and its 
intersection with Wj^ . 




Figure 1: Vi — ImA: 25000 random adapted hermitian structures 
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The modified algoritlim to estimate the shape of 7^2 — V f] A in Wj^ 
is similar. For a random R in 5*0(6), the ordered spectrum spec(C) = 
(71, ... , 76) of C = So + R^^SoR is projected orthogonally onto the diag- 
onal A by the map tta : M*^ — >■ A: 

(71,72,. ..,76) ^ (71+72^73 + 74 ^73+76)^ 

At first, when spec(C) was e-close to A i.e., if "^^t^i l72fc-i — 72fcP < 2e^ 
for e small, the projected point was plotted in green; otherwise it was 
plotted in red. No particular pattern was found for the green points 
meanwhile the red ones were all contained in Vi- The algorithm to plot 

TVAiV) is 

create 5*0 

for i = 1 to Nmax do 
create R and R^^; 
C ^ So + R-'SoR; 
1st = eigenvalues (C) ; 
sort( 1st ); 

, . f iaf[6]+isf [5] lst[4]+lstl3] lst[2]+lst[l] \ 
plot 1^ 2 , 2 , 2 J • 

end for. 

The results of the projection nA{spec{C)) for 50000 random rotation ma- 
trices are shown in Figure [2] 




Figure 2: Projection of V: 50000 random rotations 

The matrix So and hence the polytope Vi are the same as in Figure [l] 
(the interior lines correspond to the polytopes associated to different par- 
titions of the spectrum of So, as depicted in the corresponding figure in 
[CI] '). Recall that the polytope V has dimension 5; this explains that in 
order to get a better filling one should have taken many more points. This 
was not done because the evidence was sufficiently convincing to ask for 
a proof. 
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